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ADDITIONAL SOFTWARE DEVELOPMENTS WANTED FOR MODELING AND CONTROL OF 

F LEXIBLE SPACE SYSTEMS 

By 

Dr. Jiguan Gene Lin 
Control Research Corporation 
Lexington, Massachusetts 


ABSTRACT 


Existing modeling and control software packages are either Inadequate or inefficient for 
applications to flexible space structures. Some additional software developments are 
wanted for effective design and evaluation of the control systems. The following will be- .,* 
discussed, hr tilts' present aiUorT. 


1 Linear-quadratic optimal regulators as usual can be designed using various "modern 
control" design software packages. To design for active augmentation of (approximately) 

111 e sp ecified amount of active damping to each "controlled modes," the common practice is 
to adjust repeatedly the state and control weights (i.e„ the Q and R matrices) by mostly 
endless trial and error. The time consumed and effort spent in the trial-and-error 
repel it Ion can be saved by using an analytical procedure for clos ely estima ting the 
corresponding state and control weights. Various numerical examples have shewn that 
I Ids- is. possible. No software has been developed for automating such a time-saving 
analytical assignment procedure yet. 


2. "Modal dashpots" are very effective output-feedback vibration controllers for flexible 
structures, not only effective for augmenting a small amount of active damping to a large 
number of vibration modes (like the so-called low-authority structural controllers), but 
also pffeelive for Quick suppression of large vibrations (like high-authority structural 
controllers). Recent numerical results on orbital SCOLE configuration have shown so. No 
software has been developed for facilitating the design process yet. 


.1. The actual performance of any control design needs to be evaluated against a faithful 
model of the flexible structure to be controlled. The potential of destabilization or serious 
performance degradation needs to be detected by numerical simulation of the structure 
with the control loops being closed. Except for some trivial cases. reduced=order normal- 
mode models are generally not appropriate: if they are computationally feasible to 
simulate the closed-loop system, then they are likely not accurate enough to represent t e 
dynamics of the flexible structure: if they are satisfactorily accurate, then they are mostly 
loo large for effective dynamic simulation even by a state-of-the-art mainframe computer. 
Besides— eom p u ttrrg~a"veTy l a r g e numb er-of normal modes is very expensive, and the \ 

accumulated computational errors in the natural frequencies and mode shapes grow very / 
rapidly The popular Guyan reductloriHechnlque is often used to reduce the large Unite- / 
clement mass-stiffness model first. SucfLa reduction technique, unfortunately. s 
introduces large additional errors which are proportional to the square of the natural 
frequency of the modes computed thereafter. 


4 ? 


There Is a trend towards some Innovative use of non-normal modes (such as Rltz or 
Lanc7.os vectors) for representing the structures by a much smaller number of such modes 
Available results are Interesting and promising. Additional development effort Is needed 
and will be very worthwhile. 
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ADDITIONAL SOFTWARE DEVELOPMENTS URGENTLY WANTED 


• A_c curacy-Preserving Computation ally Efficient 
Cqordjln R ed u cti. on of Fi nite-Ele ment Models / 

To ENABLE 

1. Pre-Design Open-Loop Dynamic Analysis of 
Realistic, Large, Flexible Space Structures 

and 

2, Post-Design Full-Order Closed-Loop Evaluation of 
Control Systems for such Structures 


I ^^LYTjCAL_SELECnP-N o f CONTROL AND STATE WEIGHTS , 

TO AID 

Design of Linear-Quadratic Regulators desired fop 

Vibration control of Flexible Space Structures 
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ACCURACY-PRESERVING COMPUTATIONALLY EFFICIENT 
COORDINATE REDUCTION OF FINITE-ELEMENT MODELS 


Design of reliable control systems for flexible space systems 

NEEDS 

1, Careful Pre-Design Open-Loop Dynamic Analysis of 

the Space Structure, and 

2. Careful Post-Design FULL -OR PER CLOSED-LOOP Evaluation of 

Control Systems for the Structure 


Needs Pre-Design Open-Loop Dynamic Analysis 

— to assess effects of disturbances on system performance, 

E.G., POINTING STABILITY, LINE-OF-SIGHT ERRORS, ... 

-- TO IDENTIFY STRUCTURAL MODES NEEDING ACTIVE CONTROL 

— TO FORM A COMPUTATIONALLY FEASIBLE 

REDUCED-ORDER CONTROL-DESIGN MODEL 

— TO ASSSESS EFFECTIVENESS OF CONTROL ACTUATORS AND SENSORS 

Needs Post-Design FULL -ORDER CLOSED-LOOP Evaluation 

— to detect possible instability introduced by 

REDUCED-ORDER CONTROL DESIGN 

— TO VERIFY ACTUAL TIME'DOMAIN PERFORMANCE 

— TO TEST ROBUSTNESS TO MODELING ERRORS, PARAMETER VARIATIONS,,,. 
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Fig. 1-1 Spacecraft Control Laboratory Experiment (SCOLE) — 
the orbital Shuttle-Mast-Antenna configuration. 




Figure 1. A Possible Eos Polar Platform Configuration 



Figure 2. SAK Imagine fi* eroi try 
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AIAA-87-0022 paper : 

"Earth Observing System Instrument Pointing Control 

Modeling for Polar Orbiting Platform," 

by H.C. Briggs, T.Kia, S. A. McCabe, and C,E, Pell 


9 NASTRAN Structural dynamics model of Platform: 200 nodes 
-- 600 to 1,200 DOFs 
— 40 MODES BELOW 10 Hz 

I Slew of SAR antenna about its boom axis was simulated 
-- TO assess the pointing control and stability 

OF INSTRUMENTS MOUNTED ON THE CARRIER STRUCTURE 
DURING SLEWING OF ADJACENT INSTRUMENTS. 


I "Large angle time-domain simulations can presently be conducted 

USING DISCOS, BUT DUE TO EXECUTION EXPENSE AND THE DIFFICULTY OF USER 
INTERFACE THIS APPROACH IS IMPRACTICAL FOR EOS STUDIES, 

Next generation simulation tools which reduce the number of numerical 

OPERATIONS FROM ORDER N^ (DISCOS) TO N^ (TrEETOPS) AND BEYOND TO ORDER N 
ARE URGEN TLY NEEDED TO EFFICIENTLY AND COST EFFECTIVELY VERIFY 
THE PERFORMANCE OF LARGE SYSTEMS OF MULTIPLE ARTICULATED AND ROTATING 
ELEMENTS SUCH AS EOS PLATFORMS." 



SOME SERIOUS TECHNICAL PROBLEMS 


Current reduced-order models are generally not appropriate for 

REALISTIC, LARGE, FLEXIBLE SPACE STRUCTURES: 

• Model acc urac y vs comp utational. capabiliiy 

if computationally feasible to simulate on the computer, then 

LIKELY NOT ACCURATE ENOUGH TO REPRESENT THE FLEXIBLE STRUCTURE; 

IF SATISFACTORILY ACCURATE, THEN 

MOSTLY TOO LARGE FOR EFFECTIVE DYNAMIC SIMULATION ON COMPUTER 

• C 0 M P UT A I1Q NAL _E)O^NS.E_ AN D_ AC C U M U L A TED ERRORS 

Computing a very large number of normal modes is very expensive; 
Accumulated computational errors in the natural frequencies and 

MODE SHAPES GROW VERY RAPIDLY, 


9 Wasted JLXJ’fiisi ve_m od al c omp ut ations 
Many useless modes computed, 

THEN IGNORED IN CONTROL DESIGN OR EVALUATION 
— UN-RELATED TO DISTURBANCES CONCERNED, 

or Control actuations considered 

I Accuracy-S acr ificing c oor dinate re duction 

Popular Guyan reduction technique is often used first 

TO REDUCE THE LARGE FINITE“ELEMENT STIFFNESS AND MASS MATRICES 
— LARGE ERRORS INTRODUCED THEREBY; 

INCREASE AS THE SQUARE OF FREQUENCIES OR HIGHER 
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INNOVATIVE RAYLEIGH -RITZ METHOD 


I A TREND TOWARDS SOME INNOVATIVE USE OF NON“MORMAL MODES 
(SUCH AS RlTZ OR LANCZOS VECTORS) FOR REPRESENTING THE STRUCTURES BY 
A MUCH SMALLER NUMBER OF GENERALIZED COORDINATES 

-- Available results interesting and promising, 

-- Additional development and extention efforts needed, 

• Rayleigh-Ritz M ethod 

— Assumed shapes: Q2/... Q n (small n) 

-- Approximate the structural displacement vector x: 
x = z, q, + z 0 q 0 . , , + z q = 0 z 

N N 

* z ( z, , Z 0 / ,,,/ z ) 

12 N 

D D 2 x 

-- Reduce original finite-element model: M — ■ + K x = f(t) 

DT 


0 = If 


(V °2' 


TO 


(Q t M Q) ~ + (0 T K 0) z = 0 T f(t) 


Original large matrices M and K now reduced to 
smaller ones: 

M q = 0 T M 0 , K c = Q t K 0 

• Wilson- Y u a n - D i ck en s Algori thm 

— Assume f(t) = bu(t), u(t) = a scalar function 

— Generate and orthogonalize the assumed shapes sequentially: 

============> Q- 


K Q, * = B ========= 

K q * = M Ql =========> Qr 

/-I l 

« I I I • I I I I I I • I I • I I I • I I 

K Q * = M Q . ========> Q 

N N-l h 
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ADDITIONAL DEVELOPMENT AND EXTENTION EFFORTS WANTED 


| Co mputational Probl ems with Orthogona l izmion 

1, Accumulated roundoff errors can destroy the orthogonality of 

THE RlTZ VECTORS THUS GENERATED 
-- Need to re-orthogonalize whenever orthogonality is lost 

2, Computational intensive: perform Gram-Schmidt orthogonalization 

EVERY TIME A VECTOR Q * IS GEANERATED 

— Nour-Omid and Clough's solution was to orthogonalize 

ONLY WITH RESPECT TO TWO PREVIOUS VECTORS. 

— The most troublesome drawback of the Lanczos algorithm 

REAPPEAR : 

EASY LOSS OF ORTHOGONALITY OF THE LANCZOS VECTORS; 
RE-ORTHOGONALIZATION REQUIRED WHEN ORTHOGONALITY IS LOST 


I Extension beyon d th e speci al case of sc al ar forces 

— The Wilson-Yuan-Dickens algorithm was formulated for 

SCALAR FORCES; 

NOT DIRECTLY APPLICABLE TO THE GENERAL CASE OF 

MULTIPLE SIMULTANEOUS DISTURBANCE (OR CONTROL) FORCES 

— So was Nour-Omid and Clough's version using Lanczos vectors 


— But/ space systems likely be subject to multiple disturbances 

NOT ONE AT A TIME/ BUT SIMULTANEOUSLY 

— Also most control systems use multiple independent actuators 

TO APPLY FORCES/TORQUES TO THE STRUCTURES SIMULTANEOUSLY, 
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LINEAR-QUADRATIC REGULATORS (LOR) FOR FLEXIBLE SPACE STRUCTURES 


I Truncated modal model of the flexible structure 

* + 2 “ 1*1 + «2 n = ♦jBpU i = 1/ . i . /N 

I Putting into State-Space form 

x = A x + B u 


WITH 




L N J 


I LOR Design: 

Find a feedback gain matrix K such that 


J = 


'0 


(xTq x + u Tr u) dt 


IS MINIMIZED WITH U = K X 

I Given the control and state weighting matrices R and 0/ 
any "Modern Control" design program, such ORACLS, CTRL-C, 

CAN PRODUCE AN OPTIMAL SOLUTION K VIRTUALLY AUTOMATICALLY 
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DESIGN OF LINEAR-OUDRATIC REGULATORS FOR 
ACTIVE AUGMENTATION OF SPECIFIED DAMPING TO SPECIFIC MODES 


Appra oach 1 , Constrained optimization 

Optimize the performance index J with the specified damping ratios 
as constraints. 

— Constrained optimization is particularly complicated 

WHEN DYNAMIC EQUATIONS ARE INVOLVED 

Approach 2 . Alpha-shift 

Shift all poles to the left of the imaginery axis by a constant 01 . 

— Some modes may not get enough damping to be close to the specified# 

WHILE SOME OTHERS MAY GET TOO MUCH MORE THAN THE SPECIFIED. 


Appr oach 3 . Trial and error on the control and state weights 

Start with diagonal R and Q with some arbitrary numbers# e.g,# 1; 

CARRY OUT THE DESIGN OF THE CORRESPONDING LQR; 

EVALUATE THE CLOSD“LOOP POLES# AND HENCE THE DAMPING RATIOS. 

Try other control and state weights# 

REPEAT THE DESIGN-EVALUATION CYCLE# 

UNTILL THE RESULTS ARE SATISFACTORY. 

— The control and state weights used mostly are ad hoc; 

THE TRIAL-AND-ERROR PROCESS IS MOSTLY ENDLESS# 

VERY TIME CONSUMING 
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ADDITIONAL SOFTWARE DEVELOPMENT WANTED 


Software modules for aiding designers in making good initial choices, 

AND INTERMEDIATE ADJUSTMENTS, OF THE CONTROL AND STATE WEIGHTS 
SO THAT, 

THE RESULTING DESIGN OF LlNEAR-QuADRATIC REGULATORS 
CAN, WIT HIN ONLY A FEW I TERATIONS , SATISFY CLOSELY 
THE DESIGN SPECIFICATIONS, 

E.G., ON DAMPING AUGMENTATION, STIFFNESS AUGMENTATION, 

LINE-OF-SIGHT POINTING ACCURACY, ETC. 
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ADDITIONAL SOFTWARE DEVELOPMENTS URGENTLY WANTED 


I Accur acy-Preserving Co mputationally Effic ient 
Coordinate Re duction o f Finite-Element M odels , 

To ENABLE 

1. Pre-Design Open-Loop Dynamic Analysis of 
Realistic/ Large, Flexible Space Structures 

AND 

2. Post-Design Full-Order Closed-Loop Evaluation of 
Control Systems for such Structures 


I A nalytic a l Selec t ion o f Contro l and State W eights , 

TO AID 

Design of Linear-Quadratic Regulators desired for 

Vibration control of Flexible Space Structures 
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